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The effects of bulk viscosity are examined for inflationary dynamics in which dissipation and 
thermalization are present. A complete stability analysis is done for the background inflaton evo- 
lution equations, which includes both inflaton dissipation and radiation bulk viscous effects. Three 
. representative approaches of bulk viscous irreversible thermodynamics are analyzed: the Eckart 

noncausal theory, the linear and causal theory of Israel-Stewart and a more recent nonlinear and 
causal bulk viscous theory. It is found that the causal theories allow for larger bulk viscosities before 
£\J ■ encountering an instability in comparison to the noncausal Eckart theory. It is also shown that the 

causal theories tend to suppress the radiation production due to bulk viscous pressure, because 
of the presence of relaxation effects implicit in these theories. Bulk viscosity coefficients derived 

Qfrom quantum field theory are applied to warm inflation model building and an analysis is made of 
the effects to the duration of inflation. The treatment of bulk pressure would also be relevant to 
Qf) ' the reheating phase after inflation in cold inflation dynamics and during the radiation dominated 

regime, although very little work in both areas has been done; the methodology developed in this 
paper could be extended to apply to these other problems. 
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Cosmological inflationary models are usually described by the evolution of a background scalar field, the inflaton. 
These models can be separated in terms of isentropic (cold) [l|, Q and nonisentropic (warm) 0, 0| expansion (for 
■ earlier related work see [5j), when regarding the production of radiation during the inflation phase. In cold inflation 
' (CI), the couplings of the inflaton field to other field degrees of freedom are considered not to influence the inflaton 
dynamics in any appreciable way, becoming only relevant at the end of inflation, at the time of reheating. As such, 
in cold inflation dynamics, during inflation the universe undergoes a stage of supercooling. In contrast for warm 
inflation (WI) , the interactions of the inflaton field with other field degrees of freedom are strong enough to lead to 
radiation production at a rate that can eventually compensate the supercooling due to expansion. This will lead to 
ON I some general statistical state during inflation. 

The analysis of the quantum dynamics of inflationary models is done using the effective action with the in-in 
Schwinger closed-time path (or Schwinger-Keldysh) method [qj. This approach has been used for the study of the 
quantum dynamics of warm inflation since its early implementations [7|. A review of these quantum field theory 
methods applied to warm inflation can be found e.g. in Ref. 8]. In the context of cold inflation, the Schwinger- 
Keldysh method is also commonly used and one primary modification it makes to the tree level theory is to add 
time local corrections to the effective potential due to radiative effects. Another modification that can be studied 
by this method is the growth of long wavelength quantum fluctuations and their effect on the evolution of the 
inflaton background and fluctuation modes [9]. Warm inflation models are those in which significant time nonlocal 
contributions also arise from the effective action. Main consequences of such terms will be dissipation and fluctuation 
effects @, • This can very generally lead to a vast variety of particle production behavior during inflation, thus 
to various types of statistical states. Amongst all of these, the one case which can be studied most extensively by 
quantum field theory methods is the case where the statistical state of the Universe emerging from these effects is 
near thermal equilibrium. Once particle production occurs during inflation, there are many effects that can follow 
and it is a complicated problem. This is why it is essential to understand the properties of warm inflation in the near 
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thermal equilibrium regime, so as to provide insight into the features that one should also expect when away from 
this limit. Of late there have been some efforts to understand warm inflation models away from thermal equilibrium 
fill . These papers have analyzed some of the effects that emerge in the presence of particle production, but they 
are at a much more primary level then the extensive study for the near thermal equilibrium regime. In particular, 
creating a radiation bath of particles will also then modify the inflaton potential no matter what statistical state it is 
in. Also the nature of the radiation bath of particles will alter and modify the density perturbations. In the thermal 
limit there are well established methods for examining all these effects, whereas away from this limit it is a much 
more difficult problem. 

As mentioned above, the case that has been studied most extensively, since it is the most tractable by quantum 
field theory methods, is when the radiation has adequate time to reach close to a thermal state. For this case, the 
radiation production during inflation manifests itself through an effective dissipation term in the inflaton background 
equation. The corresponding dissipative coefficient can be computed from first principles in quantum field theory 
under appropriate conditions (see e.g. Refs. IS GIHIII] ) ■ In particular, the two-stage interaction configuration for 
quantum fields proposed in [3, [13], has been shown can lead to a large enough dissipative terms, while ke eping t he 
quantum corrections to the inflationary potential under control and allowing a period of slow-roll inflation [18l - l20l | . 

In WI, the transfer of energy of the inflaton field to the radiation bath is mediated by the dissipation term in the 
inflaton's evolution equation. Nevertheless, an additional effect can arise due to inner couplings in the radiation fluid 
itself. Internal decays within the radiation fluid can make it depart slightly from thermal equilibrium. Therefore, the 
radiation fluid can behave as a non-ideal fluid and viscosity effects must be taken into account [2l| . At the background 
level, the relevant viscous effect is due to bulk pressure, since it is the only viscous effect appearing in the background 
equations for an FRW universe. 

The study of bulk viscous effects in cosmology, and in particular in inflation, has some history to it, focused mainly 
on the effect of the bulk pressure as a negative pressure (for a partial sample of the earlier works on bulk viscous 
cosmologies, see for example Ref. [22T]). In addition, more recently, there has been a surge of interest in exploring 
the effects of the bulk pressure as the origin of the present accelerated expansion of the universe (see e.g. Refs. [Hj]). 
Almost all of these past works only use phenomenological forms for the bulk viscosity. This paper differs from these 
past works since we will apply first pri ncip le quantum field theory computed expressions for the bulk viscosity, based 
on the calculations found e.g. in Ref. [24| . to warm inflation. There has been very little work done in studying first 
principles bulk viscosity expressions in application to cosmology. There is one paper we are aware of along these 
lines, which is the first reference in [23|; this paper uses quantum field theory derived expressions of bulk viscosity 
and shows how they can play the role similar to dark energy. 

Our work will examine the consequences of bulk viscosity on warm inflationary dynamics. In particular we are 
interested in determining the stability conditions of the background equations of warm inflation when coupled to 
the bulk viscous radiation bath. Earlier studies examining the stability properties of the warm inflation equations 
include Refs. |20L l25j. but these studies did not include the effects of bulk viscous pressure. Preliminary studies on 
the inclusion of bulk viscous pressure have been done in [26[ , [27j , where only the non-causal theory of Eckart [28| has 
been used. Here we extend the stability analysis of the dynamical warm inflation e quat ions to include bulk viscous 
effects. Moreover our analysis will be done not just for the non-causal Eckart theory [28[ but also the causal theories 
of bulk pressure [29M3ll j . with a full analysis of the differences in the resulting dynamics from these different theories. 

In the study of viscous effects in cosmology it is common to use linear expressions to describe the viscous pressure, 
where it is assumed that the deviations are close to equilibrium. However, the viscous pressure can, in principle, 
take the system far from its thermodynamical equilibrium, so, we must apply suitable approaches in order to see if 
one really needs to use a more robust description, incorporating nonlinear effects. We consider here three different 
theories to describe the viscous pressure: the non-causal theory due to Eckart [28], the linear and causal theory of 
Israel-Stewart (29l . l30j and finally we will also use a recent causal and nonlinear theory proposed by the authors in 
named by them Nonlinear Causal Dissipative Hydrodynamics (NLCDH). There have been other proposals for a 
nonlinear theory for the bulk viscous pressure [H, l3p. that make use of ad-hoc parameters, such as the time where 
nonlinear effects become important [32J or functions [33] that do not have an immediate interpretation from quantum 
field theory, for example being a relaxation time. In using such approaches there is no immediate understanding how 
to associate their parameters with first principles parameters. We have considered the theory for bulk pressure in [3l[ 
since it utilizes parameters which can readily be determined from microscopic physics, in particular from quantum 
field theory. The theories for the bulk pressure we analyze here are more natural to use in field theory model building, 
where the dissipation terms, viscosity coefficients and relaxation times are well defined and can be reliably computed 
once a specific field theory model is given. 

We study the effects of the inclusion of bulk viscosity in three commonly used supersymmetric realizations of warm 
inflation, the chaotic, hybrid and hilltop models. As the bulk viscosity modifies the background dynamics of warm 
inflation, it also changes the available parameter space for warm inflation, which is analyzed here. The bulk viscosity 
is approximated by the Eckart description and we place limits on the couplings of the underlying particle physics 
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theory for the validity of this approximation. 

The paper is organized as follows. In Sec. [TT] we summarize the relevant expressions for warm inflation and for 
the theories of bulk viscous pressure that is relevant for this work. In Sec. Mil we present the respective dynamical 
systems in each of the three theories for the bulk viscous pressure we are examining, along with the resulting stability 
equations. In Sec. IIV1 we present some numerical results that confirm the stability results. Also in this section there 
is an exposition about the differences in the dynamics that emerge when considering noncausal and causal theories for 
the bulk pressure. In Sec. [V]we discuss quantum field theory derived expressions for the bulk viscosity and conditions 
under which they can be applied to cosmology. In Sec. IVII we apply the stability analyses to explore the parameter 
space available for warm inflation constructed from supersymmetric models. Finally, in Sec. IVIII we present our 
conclusions and final comments. 



II. WARM INFLATION IN A BULK VISCOUS RADIATION FLUID 

In the warm inflation scenario, the inflaton evolution is described in terms of a dissipative equation of motion of 
the form 0, 0,0 El 

0+ {3H + T)<j> + Vj, = 0, (2.1) 

where T is the dissipative coefficient, is the field derivative of the (field and temperature dependent) effective 
potential, H = a/a is the Hubble rate of expansion, and a the scale factor of the Friedmann- Robertson- Walker 
background metric, ds 2 = —dt 2 + a(t) 2 5ijdx l dx : > . 

The dissipation term in Eq. (|2.1j) originates from the interactions of the inflaton field with other fields. These 
interactions effectively render the inflaton as an open system, which necessarily should give origin to dissipation. 
The dissipation coefficient T accounts for this effective dissipation of the inflaton and, given a particle physics model 
realization of inflation, it can be computed microscopically (see e.g. the recent review Q and references there in). 
If enough dissipation, i.e., radiation, can be produced during the inflationary regime, it can eventually compensate 
the redshift of the radiation due to expansion. When this happens, we then term the inflationary regime as warm 
inflation. A particular field theory model that has been proved successful in generating the WI conditions is the 
two-stage mechanism [l6l - fl8l |. In the two-stage mechanism, the inflaton field couples to a heavy catalyst field, and the 
latter in turn couples to light degrees of freedom. The evolution of the background inflaton field excites the catalyst 
fields, which then decay into light fields. Large dissipation has been demonstrated with this simple mechanism, 
adequate to drive a warm inflationary expansion. When the produced light relativistic particles thermalize, which is 
achieved when their creation and scattering times are less than the Hubble time in an expanding universe, we can 
model their contribution as that of radiation, 

p r ~ ^.g,T 4 , (2.2) 

where T is the temperature of the thermal bath, and <?» the effective number of light degrees of freedom 1 . 

Bulk viscous effects come in to play when the light particles in the resulting radiation bath in warm inflation can 
also decay, thus making the radiation fluid slightly depart from equilibrium. This effectively renders the radiation as 
a nonideal fluid and gives origin to a bulk viscous pressure (there may be other dissipative effects in the radiation 
fluid itself, like shear viscous stresses, but these do not contribute at the background level [2l| . though they may 
be relevant in the determination of the spectrum of density perturbations, as demonstrated recently p^V In the 
presence of a bulk viscous pressure, n, the stress-energy tensor for the radiation fluid is given by [2ll l35j|. 

rf2 = (Pr +Pr+ n)«W«W + ( Pr + H)^ , (2.3) 

(r) 

where p r is the radiation energy density, p r the adiabatic radiation pressure, are the four velocity of the radiation 
fluid and g^, the four-dimensional metric. It happens then that the bulk pressure enters as a contribution to the 
radiation pressure p r , such that we can define in general an effective pressure for the radiation, p r , given by 



If not otherwise specified, we will take g* = 228.75, the effective no. of degrees of freedom for the Minimal Supersymmetric Standard 
Model, when presenting numerical results. 
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p r = p r + H. (2.4) 
The evolution equation for the radiation fluid energy density p r then becomes 

p r + 3H(p r +p r ) = T(p (t> + p<j,) , (2.5) 

where p$ = <j) 2 /2 ~ V(4>, T), and pj, + p^ = 4> 2 . 

It is also useful to express this in terms of the entropy density s. From the Helmholtz free energy f — p — Ts, where 
/ = V(cj),T), and using s = —df/dT, the total energy density p becomes 



P=^+V(4>,T)+Ts, (2.6) 



and the Hubble rate H becomes 



H 2 = 



3m| 



^-+V(cf>,T) + Ts 



(2.7) 



where mp is the reduced Planck mass, rap = l/v87rG = 2.4 x 10 18 GeV. Using also that p r = (7 — l)p r and that the 
entropy density s is related to the radiation energy density by Ts = 7p r , Eq. (|2.5[) can then be written in terms of 
the entropy density as 



Ts + 3H(Ts + 11) = T<j> 2 , (2.8) 

where we have used 7 = 4/3, which is valid for a quasi-equilibrium high temperature thermal bath typical of warm 
inflation. 

From Eq. (|2.8[) . the dynamical effects of the bulk viscosity can be easily read. Given that the bulk viscous pressure 
II is negative, it acts to decrease the radiation pressure, thus enhancing the effect from the source term on the RHS 
in the equation for the entropy density. As a consequence, the entropy density increases, and therefore the radiation 
energy density also grows. On the one hand, if this bulk pressure term is too large, there is too much radiation 
production and the radiation energy density dominates too soon over the scalar field energy density, thus spoiling 
inflation. This regime is called the unstable regime. On the other hand, if the bulk pressure term is controlled to 
avoid the radiation domination until the end of inflation, the system is said to be in the stable regime. In this regime 
the bulk viscosity gives rise to an additional negative pressure, and hence, inflation is enhanced. 

To account for the dynamics involving the bulk viscous pressure II, as explained in the introduction, we will consider 
three different theories: the noncausal theory due to Eckart [28}, the linear and causal theory of Israel- Stewart 
(IS) (23 . |30|, a nd a recent causal and nonlinear theory, Nonlinear Causal Dissipative Hydrodynamics (NLCDH), 
proposed in [3l[ . The starting point to build these hydrodynamic theories is the conservation equations of the stress- 
energy tensor and the number density vector N 11 = nu^, 

V M T^=0, V fi N* = 0, (2.9) 

with the additional condition on the 4-entropy, written in terms of the entropy density s, s M = su^, that must satisfy 
the second law of thermodynamics in its covariant form, 

V^>0. (2.10) 

The 4-entropy, just like the stress-energy tensor, gains a contribution coming from the dissipative fluxes, 



fl f = « t " + ^_ i (2.11) 

where Q^ 1 — Q^(N^, T MI/ ) accounts for the dissipative fluxes. The irreversible thermodynamics comprises of the 
dissipative forces to the hydrodynamics variables at equilibrium, the number density n, the energy density p and the 
pressure p. These quantities are able to describe the energy fluxes in a nonideal fluid. There are different ways in 
which this can be done, which lead to different descriptions for the dissipative fluxes, like for example for the bulk 
pressure. We summarize below the Eckart, IS and NLCDH theories for the bulk pressure. 
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A. Eckart theory for the bulk pressure 



The Eckart theory 28| assumes that the entropy vector s 11 is linear in the dissipative fluxes. The nonequilibrium 
contribution to the entropy vector, Q M , to first order, should then be proportional to the dissipative fluxes. Neglecting 
dissipative terms other than the bulk pressure, we then have that 



Q* 1 ~ a(n, pjHuf* , (2.12) 

where the proportionality factor is obtained from the equilibrium condition and from the covariant form of the second 
law of thermodynamics, Eq. (I2.10j) . This then gives [35| 



TV^ ~ -3HU . (2.13) 

To ensure that the second law of thermodynamics, Eq. (|2.10p . is satisfied and interpreting the term 3H in Eq. (j2.13p 
as a dissipative force, \e = 3H, we impose II to be linear in this force, expressing bulk viscosity as 



Il = -3C b H, (2.14) 

where the proportionality term = p) > is the bulk viscosity coefficient j2ll |. 

The bulk pressure expressed like Eq. (|2.14[) is a noncausal theory, i.e., the speed of the fluxes propagation is 
infinite. The Eckart theory can be considered in some circumstances as a reasonable approximation for the irreversible 
thermodynamics. This may happen, for example, when sufficiently short relaxation time scales are involved, otherwise 
a causal theory would be a much better choice. We now turn to the simplest of such a causal theory, the IS one. 



B. Israel- Stewart theory for the bulk pressure 



The IS theory [29|, [30( goes one step further than the Eckart theory by accounting for second order contributions 
beyond equilibrium, by expanding the entropy vector to second order in the dissipative fluxes. Generically this gives, 
by again only considering the bulk pressure contribution, 



S ^_/3 II 2 — , (2.15) 



where /3a(n,p) > 0. From the covariant derivative of the entropy vector, 



-n 



3H + /3 n + ~v„ (^u» ] n 



(2.16) 



and from the second law of thermodynamics to be satisfied, Eq. (|2.10l) . it is imposed again, like in the Eckart case, 
that the dissipative fluxes be linear in the dissipative forces. For the bulk pressure II this implies from Eq. (|2.16[) that 



(2.17) 

The relation (I2.17[) is analogous to Eq. (|2.14l) in the Eckart theory. The difference here being that, from Eq. (|2.13p . 
the Israel-Stewart force is expressed as 



n = -c& 



3H + M + |v„ f |y l n 



Xis = + ATI + |V„ n . (2.18) 

By defining r = CfcAh which is interpreted as a relaxation time for the bulk viscous processes in the radiation fluid, 
then Eq. (|2.17p can also be rewritten in the form 
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and by expanding the derivative in the last term in Eq. (|2.19p it can finally be expressed as 

T n + n = -3C6ff-^^ + i-|-^ , (2.20) 

which is the IS equation for the bulk pressure. 

As shown in Ref . [3(| , the propagation speed for the bulk pressure is given by 



c visc — / . \ ) (2-21) 

(P + PF 

and, thus, for t / there is a finite propagation speed for the flux, while for the Eckart theory, where r = 0, it is 
infinity (noncausal). In a quantum field theory description for the radiation bath, e.g. for example in the two-stage 
decay mechanism mentioned earlier for warm inflation, both the bulk viscosity coefficient and the relaxation time r 
can be defined unambiguously and be computed microscopically, just like the dissipation coefficient T. In particular, 
the bulk viscosity coefficient can be obtained from a Kubo formula for the high-temperature light particles of the 
radiation bath [24], and r can be associated with the respective decay time of these particles, t = 1/r, where T is 
the decay width. For the validity of considering a quasi-equilibrium thermal radiation bath, we are then required to 
impose that 



tH = H/T < 1 . (2.22) 

Likewise, the assumption of proximity with thermal equilibrium requires the dissipative fluxes to be small compared 
to the equilibrium pressure, 



|n|<p. (2.23) 

In terms of the quasi-equilibrium conditions (I2.22[) and (I2.23[) . the equation for the bulk pressure can also be 
expressed in terms of the Maxwell-Cattaneo equation for II |37| . 

rri + n = -3( b H. (2.24) 

The IS equation for II, Eq. (I2.20p . can then be seen to give a correction to the Maxwell-Cattaneo equation, which, 
after Eckart, are the simplest equations for the bulk pressure including relaxation (causal) effects. 



C. Nonlinear causal dissipative hydrodynamics theory for the bulk pressure 

Next, let us consider the NLCDH theory proposed by the authors of Ref. [3l|. This theory assumes the Eckart 
force term, \e — V M u M = 3H , plus a memory effect, so as to respect causality. Since this theory in principle makes 
no assumptions about the linearity of the bulk pressure, as it was assumed in the IS theory for instance, it has been 
regarded as a nonlinear theory for the bulk pressure. The memory effect adds a relaxation to the system. Recall 
that the Maxwell-Cattaneo theory Eq. (|2.24[) is obtained by adding a relaxation time directly to II. In the NLCDH 
instead, the memory effect is added to the quantity II = LIV (where V is the volume), which is then integrated in a 
cell (volume element) of the fluid. This is done by imposing the relation II = — Cb^XE- Through the addition of the 
memory effect, we are lead to [3l[ 



tH + 11= -3H( b V. (2.25) 

The first term in the above equation results in r(IIV + ITV). After using the conservation law for the volume in a cell 
of the fluid, V M (u^/V) = 0, it can be shown that [U 
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v*Y = {v) + v v ^ = -{v) v+ v* = °- (2 - 26) 

It follows that x = V/V and then V = %V = 3HV. From Eq. (|2~25j) . it then follows the NLCDH equation for the bulk 
pressure [3l| . 

tII + n = -3H(Cb + tIL) . (2.27) 

III. THE DYNAMICAL SYSTEM OF EQUATIONS FOR WARM INFLATION IN A BULK VISCOUS 

RADIATION FLUID 

The relevant equations concerning warm inflation in a bulk viscous radiation fluid are given by the inflaton evolution 
equation, Eq. (|2.1I) . the entropy energy density evolution Eq. (I2.8p . with the bulk pressure II given by: (a) in the 
Eckart case Eq. (|2. 14|) : (b) in the IS case by the evolution equation ()2.20j) : and (c) in the NLCDH case by the evolution 
equation (j!T2"7]) . 

Writing the inflaton equation of motion as two first order differential equations, we have that Eqs. (|2.ip and (|2.8[) 
are equivalently written in the form: 

</» = u, 

ii = SHu — Tu—V^, 
Ts = -3HTs-3HU + Tu 2 , (3.1) 

where 



3^{Y + V + Ts ) 1 V = V ^ T ^ T = Tfo,T), (3.2) 

with expression for the bulk pressure II, given by either Eq. (|2.14l) , (|2.20p or (|2.27p , depending on which of the cases 
is treated. In all the cases, they depend on the bulk viscosity coefficient, Cb = 0>(T). In warm inflation the expected 
dependence of the dissipation coefficient on the field and temperature is [si ll3l.Tl5| 

r = C^, (3.3) 

with proportionality factor depending on the field content of the model and the value of the power c depending on 
the temperature regime for the different fields involved. For example, in the two-stage decay model for warm inflation, 
for a radiation bath of light particles y in the high temperature regime, m y <C T, and for heavy catalyst fields x m 
the low temperature regime, m x > T, we find that c = 3 fl3l [20|. Likewise, the bulk viscosity coefficient, in this 
same regime is [2~ij £5 oc T 3 . In the following we consider a generic power dependence I in the temperature for the 
bulk viscosity coefficient, C,b oc T, similar to that considered for the dissipation coefficient, Eq. (|3.3|) . 

Treating the variables of the dynamical system in the form of a column matrix X, we can express the dynamical 
system in the compact matrix form, 

X = ¥{x) X , (3.4) 

where for example, for the dynamical system given by Eq. (|3.ip . X = (cf>, u, s). In the IS and NLCDH cases we also 
have the bulk pressure entering in the system as an additional function, X = (0, u, s, n). 

Writing x = xq + 5x, where we assume that xq is a stable solution of the system, which here will be taken as 
the slow-roll solutions that can be derived directly from Eq. (|3.1I) (see below), the equation for the variations in 5x 
become 

S± = M(x ) SX - X , (3.5) 



8 



where 

M( X0 ) = , (3.6) 

is the Jacobian matrix for the system, evaluated at the Xq solution, and X is a residual force term, which in general 
is small and can be neglected |2(j . The general solution of Eq. p. 51) is of the form 

<5X = Xoe^ ' 4 , (3.7) 

and M(iEo) must be zero or negative for the system be stable, i.e., the eigenvalues A, of M(xo) must necessarily all 
satisfy 

At < . (3.8) 

The stability of the dynamical system can be studied directly in the time variable, but it simplifies the analysis, in 
particular the determination of the eigenvalues of the Jacobian matrix, if we make a change of variables and rewrite 
the dynamical system (|3.1[i using <f> as the independent variable instead of the time [20| . By doing this the dimension 
of the corresponding system is smaller and easier to analyze. In particular, the Jacobian matrix for the dynamical 
system (|3.1[) becomes a 2 x 2 matrix; if one includes the bulk pressure as an additional function to the system, as in 
IS and NLCDH cases, it then makes the Jacobian matrix 3x3. The eigenvalues obtained by using <p rather than time 
as the variable makes the analysis much simpler. As such, using that 

d_ = f#_£ = tt -l = u ()' (3 9) 

dt dt dip dtp 

where a prime indicates derivative with respect to <p, the dynamical system Eq. (|3.ip becomes equivalent to 

u = -3H - T - V^u' 1 , 
Ts' = -3HT.su- 1 - SHUu- 1 + Tu , (3.10) 

together with the corresponding equations for the bulk viscosity, Eqs. (|2.14[) . (|2.20p or (|2.27|) (these last two also 
transformed to the <p variable), depending on which case we are considering. 
The general solution of Eq. (|3.10p is now of the form 



SX = X e M(KoWt) . (3.11) 

The stability condition on the eigenvalues Ai, which are the eigenvalues of M(xo) once u is factorized, depends now on 
whether the inflaton field <p(t) during slow-roll decreases with time (like in chaotic inflation), so A^ > 0, or increases 
with time (like in hilltop inflation), in which case Ai < 0. This is an important consideration when replacing the 
time by the inflaton as the independent variable in the dynamical system. In a 2 x 2 system, stability is ensured 
once the determinant is positive and the trace, negative, as is done in [2(|, (26|. However, in a 3 x 3 system, the case 
for the Israel-Stewart and NLCDH descriptions, further information beyond the trace and determinant is required 
to guarantee the stability of the system. That is why in this work we will make use directly of the eigenvalues for 
searching for the stability of the system. Nevertheless, for all the cases we have studied below, actually only one 
eigenvalue changes sign at the instability point, and therefore it would be enough to look at the determinant of the 
system. 

As mentioned above, the stability analysis is done around the slow-roll solutions xq. We study the stability of 
dynamical svstem around the slow-roll solutions, since they act as formal attractor solutions for the dynamical 
system [25|, [38[ . The slow-roll solutions are defined by the equations 

V,4, 

U 3H(1 + Q) ' 
Ts = Qu 2 - n , 

n = -3( b H , (3.12) 
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where we have omitted the index that indicates the slow-roll solutions used in Eq. (13.51) for simplicity. Note that 
we are writing the slow-roll condition for II as the same in all three cases for the bulk viscous pressure that we have 
considered in this work. This is valid provided that the relaxation time r is sufficiently small, which can be implicitly 
assumed by the condition (|2.22p . In Eq. (I3.12[) we defined Q = Y/(3iT), with the Hubble rate H, in the slow-roll 
approximation as given by 



H 2 = — t V(1 + k), (3.13) 
3mp 

where k = p r /V. Keeping the radiation energy density in Eq. (|3.13[) is justified because in the presence of a bulk 
pressure, the radiation energy density does not in general need to be much smaller than the vacuum energy density 
in order to have inflation. In particular recall that in the presence of a bulk pressure, the acceleration equation is 



6mp 



(2V + 3|n| - 2p r ) , (3.14) 



which shows that we could in principle have p r ~ V and inflation could still be sustained by the bulk pressure 133 
The solutions (|3.12j) follow if the following conditions on slow- roll coefficients in warm inflation are satisfied [20| , 



p 

~ V v 



v - 2 



V 



«1 + Q, 
< 1 + 0, 



where the slow-roll parameter f3 takes into account the variation of T with respect to cf>, and the last condition for b 
ensures that thermal corrections to the inflation potential are negligible. 

Let us now study the stability of the dynamical system of equations in each of the three cases for the bulk pressure 
that we are considering. 



A. Dynamical system for the Eckart case 

In the Eckart theory the bulk pressure is simply given by Eq. (|2.14|) . II = — 3£&-ff. The dynamical system Eq. (|3. 10|) 
takes the form, 

u' = -3H -T -V^u- 1 = f(u, s) , 

s' = -3HSV,- 1 + 9( b H 2 (Tuy 1 + TT~ 1 u = g{u, s) . (3.16) 
The Jacobian matrix M becomes 



UEckart 



(X) 



d(f,g) 



d(u, s) 



= / df/du df/ds 
~ \ dg/du dg/ds 



U—Uq,S — S() 



U — Uq ,S — Sq 



A S 

e d 



(3.17) 



where the matrix elements are evaluated at the slow-roll solutions Eqs. (j3. 12[) and (|3.13[) . The coefficients of the 
matrix Msckart become 
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A 



D 



H 

s 



3(7 -l)oQ + 3(7 -1)6(1 + Q) + 



1 



(l + /t) 2 (1 + Q) 2 1 + k 
1 e 



1 



U 

H 

u 



6 



(1 + k) 2 (1 + Q)- 



(1 + s) (1 + Q) 

f [6(l + Q) 2 -2 £ ] 
< 1 + cr- 



3( 7 -l)(c-l)-3- 1 ^ 



[6(1 + Q) 2 - e] 
Qe 

(1 + K ) 2 (1 + Q) 2 



3( 7 -l)(c-0- 



1 



Qe 



3 cr 



(1 + k) 2 (1 + Q) 2 2 1 + k_ 
where we have omitted the sub-index "0" of the slow-roll solutions and defined the quantities a and cr as 



(3. IE 



n 

IPr 

n 



(3.19) 
(3.20) 



The expressions simplify considerably in the strong dissipative regime of warm inflation, Q > 1 and neglecting the 
terms proportional to the slow- roll parameters in Eq. f|3 . 18[) , we obtain 



H 



A = -3Q- 

a 

5 = 3( 7 -l)(6-c)Q 

6 = ^ 6(1 + a), 



H 



V = -^3( 7 -l)(c-l)-3- 



3(7-D(c-0 + ^ 



Using p.21[) . the eigenvalues of M Eckart are 



Eckart 



- [3Q + 6(1 + a)(7 - 1)(6 - c)} + (1/Q) , 



v Eckart 



H 



3( 7 -l)(c-l)-3 + 



3( 7 -l)(c-0 + 



21+K 



+ 6(l + £r)(7-l)(6-c) } + 0(1/0) 



Independent of the inflaton dynamics, we then obtain that stability is assured when (u/ H)\ < 0. 
Note that in the slow-roll regime we have for cr = 11/ ("fp r ) that 



(3.21) 



(3.22) 



(3.23) 



n 



Qu 2 -n 



2Q e__ 

(1 + k)(1+Q) 1+Q 



(3.24) 



where we have used the slow- roll equations for the radiation energy density and u, Eq. (|3.12p . Note from the above 
equation that in particular we have that \a\ < 1. Using (|3 . 24[) in Eqs. (|3.22l) and p.23| . the first eigenvalue will 
always satisfy the stability condition, while for the second eigenvalue Eq. ()3.23j) . the stability condition implies: 
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(c - 26) (1 + „) + ^ + la - > , (3.25) 

or, using 7 = 4/3, valid for the quasi-equilibrium thermal bath of warm inflation, 

(c-2h)(l + <7) + 4 + Z<7- > . (3.26) 

2 1 + K 

The Eq. (|3.26|) generalizes the results in 26], which were obtained for a constant bulk pressure (1 = 0), by accounting 
for temperature dependence 2 . Also for a = and a = 0, the case of zero bulk pressure, we reproduce the results 
obtained by Moss and Xiong in [2(|. In [2(| the stability condition was found to be |c| < 4. From Eq. (|3.26j) . in 
the absence of bulk viscosity, we derive instead only the condition c > —4. We do also obtain the result c < 4 if we 
consider the eigenvalues in the approximation of very small dissipation Q <C 1, but this regime is not the most general 
situation for warm inflation |39| . 



B. Dynamical system for the IS case 

Let us now consider the dynamical system when the bulk pressure II has an evolution according to the IS theory, 
Eq. p.20p . The dynamical system, including the evolution equation for the bulk pressure, now becomes 



-3^-T-y^- 1 = f(u,s,IL) 
SHsu- 1 - SHU (Tu)~ 



n' = 



n 

u 

T 



3( b H _j n I 

u < itlu 

2 I 



T 



r 



- (7-1 

and the Jacobian stability matrix becomes 



Cfc.T 



g(u, s, II) 



t Ct V T Cb 



V 



3HU- 1 1 



s 

n 

Ts 



Tit 
7^" 



h(u, s, n) 



(3.27) 



M IS (x) = 



d(f,g,h) 



d(u, s, n) 




(3.28) 



u—uq,s—sq ,ri— rio 

Using the slow-roll solutions, Eqs. (|3.12|) and (|3.13j) . we obtain for the elements of the matrix M/s the result 



2 Our result given by Eq. 13.2611 also corrects Eq. (22) in 26] , where the authors of that reference mistakenly replaced a by a. 



12 



A 



H 

u 

H 

s 



- 3(1 + 0) -(TW(TTqF. 

- 3( 7 - l)cQ + 3(7 - 1)6(1 + 0) + 



1 



= E 

u 

£ = , 



(1 + k) 2 (1 + Q) 2 ' 1 + K 

1 e 



1 



(1 + k) (1 + Q) 2 2* 



6- 



(1 + a) , 



(1 + k) 2 (l + 0) 2 _ 
3(7-l)(l + ff)c-3-3(7-l) 



1 



(1 + K ) 2 {l + Qf 



'J 



H 

> 

Tu 



S = -^2- <7 { 2 + 3 (7-l)A(l + <r)+ 



+ 



1 



3(1 + nf 



1 3 3(7-1) 
9 2 2 



A(l + a) 



ffT f ( 7 -l)/ , 3( 7 -1) ; 

(7 



+ 



// 1 e 

3(7-1), (1 + 0) 



3(1 + /c) 





1 3 

9^2 



2 

(1 + *) 



cA(l + a) 





(1 + O) 2 

37(7 - 1) 



A- 



1 + 

3(7- 1)A(1 +a) 



A + (7- 1)E 
1 



(7-1) X 

2(1 + «) 1 + Q 



Oc 



3. 



3= ^ 

it 



1 3 3(7-1) 



e 2 

where we have defined the parameters 



A 



1 + k(1 + 0) 2 2 

1 3-7 



2(1 + /s) 1 + 



(3.29) 



Q = tH , 1= T ^ T , 1= 1 {Cb ^ 

IV T - Ttst _ 

.9 = — - , 9 = ) 9 



8nG V CbV 

1 / T 



b,TT 



9 



(b,T 

Tt tt 



t T0 8nG \ tV J t,t 

A = 1 + 1 - g , x = (g-g)g-{i-i)h 

E = (l+g- g)g- (1 + 1-1)1- Ab, b= TV>rA> • 



(3.30) 



In the strong dissipation regime 0^1 an d neglecting terms proportional to the slow-roll coefficients, Eq. (|3.29j) 
simplifies considerably and we obtain 
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A 
B 

e 

D 
£ 

s 



-3Q— , 

u 

3( 7 -l)(6-c)Q- 
s 

^6(1+<t), 



u 
0, 



[3( 7 -l)(l + a)c-3-3(7-l)] , 



H 

HTs 



u- 
HT 



+ 3(7 - 1)A(1 + a) 



(7-1)1 , 3( 7 -l) 1 



3(7-1) 



cA(l + o-) 



37(7-1) 



2 



2(1 + k) 

u \ e 2 

In terms of Eq. (|3.32[) , the eigenvalues of M75 are 



6(1 + 0-) [A+( 7 -l)S] + 
I-|-3( 7 -l)A(l + «7) 



3(7-1) 



A [a + 6(1 + a)] 



A- 



(3.31) 



(3.32) 



is 



A 18 

A 2 



A 3 



-- [3Q + 6(1 + ff )( 7 - 1)(6 - c)] + (1/Q) , 











1 






») 










+ 2 






se 




1 






») 








"X " 


~ 2 



f'B - :) ! ( D — J — 2^) - i:r f ^ - jc 



(D-5) D-J-2 



~A~ 



+T I -Jl_ J£ 



s 2 e 2 

■s 2 e 2 
Ta 2- 



1/2 



1/2 



(3.33) 

+ 0(1/Q), (3.34) 
0(1/0) • (3.35) 



For the first eigenvalue above, A', we see that it is the same as the Eckart case, Af ckart , so it already satisfies the 
stability requirement that (u/H)Xi < 0. Adding all the factors in the two other eigenvalues in Eqs. (|3.34[) and (|3.35|) . 
using that 8 <§C 1, we obtain that 



A 



is 



6 



3H 



(7-1)' + 



2(1 + *) 



0(6) 



(3.36) 



and it also satisfies the stability requirement. The stability condition for the IS case then falls on the third eigenvalue 
A| s , Eq. (13.351) . This is most easily expressed by demanding that the product AJfAg 8 > 0, which then leads to the 
condition 



2cr + 3(7-l)e6A(l + cr) 2 



2 + 39 



c + 



7 



2a 



-I 



7-1 2 + 36 



00 



2(7-1) 1 + K 



2-36 6( 7 -l)6A 



2 + 36 



2 + 36 



(1 + t) 



4(1 + a) + 36(2 + a) - 36A(1 + af + 

3(7 - 1)6(1 + cr)(A + ffS) + 66(7 - 1)A6(1 + cr) 2 



2 + 36 



>0, 



(3.37) 



which for 7 = 



4/3 becomes 
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2a + 06A(1 + a) 2 



2 + 39 



2a 



36 



Z 3 aa 



21 + K 



2-39 



29A 



2 + 39 2 + 39 



(1 + *) 



4(1 + a) + 39(2 + cr) - 39A(1 + a) 2 - 9(1 + a) (A + <r£) + 2feA9(l + a) 



2 + 39 



> . 



(3.38) 



From Eq. (|3.38[) , when the relaxation time vanishes, r = 0, i.e., for 9 = tH — > 0, we recover the previous condition 
Eq. (|3.26[) , obtained in the Eckart theory case. 



C. Dynamical system for the NLCDH case 



Finally, we will now obtain stability condition for the case of the NLCDH theory for the bulk pressure. In the 
NLCDH case, the evolution equation for the bulk pressure is given by Eq. (|2.27p . The dynamical system now 
becomes 



u = -ZH -T -V^u' 1 = f(u,s,U) , 

s' = -SHsu- 1 -SHUiTuy 1 + TT- 1 u = g(u,s,U) , 

n' = --u- 1 -^iu- 1 -SHUu- 1 = h(u,s,Il) . (3.39) 

T T 

The Jacobian stability matrix is similar to the one in the IS case, Eq. ()3.28j) . but now with the functions 

f(u, s,n), g(u, s,n), h(u, s,II) obtained from the above equation (|3.39[) . Using again the slow-roll solutions, 
Eqs. p,12p and (|3.13l) . we obtain for the elements of the matrix Mnlcdh in the NLCDH case, 



A 
■B 



3 



-3(1 + Q) - 



(i + K ) 2 (l + Q) 2 



H 

u 

j | - 3( 7 - l)cQ + 3( 7 - 1)6(1 + Q) 



1 



(1 + k) 2 (1 + Q) 2 1 + k 



Hs 

„,2 



6 



(1 + k) 2 (l + Q) 2 







T> = 






u I 


£ = 


0, 


J = 












HTs 


S = 


2 CT 



.(1 + k) (l + Q) 2 
(1 + » 

1 



3 _' 
—a 
2 



(l + /t) 2 (l + Q) 2 



i?T f (7 - l)i 
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1 / 1 3 

3(1 + k) 2 ^9 ~ 2; (l + Q) 2 
1 3 



1 



9 



-9- 3 



3(1 + k) \e 2 



Qe 



(1 + k) (l + Q) 2 



3_ 
—a 
2 



(3.40) 



with parameters the same as defined in Eq. (|3.30l) . 

Taking again the strong dissipative regime, Q > 1, and neglecting terms proportional to the slow-roll coefficients, 
Eq. (I3.40|) simplifies to 
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A 
S 

e 

D 
£ 
7 

S 
"K 

3 



-3Q- , 

u 

3( 7 -l)(&-c)Q- , 
s 

Hs 



ti- 
ll 



6(1 + a) , 



[3( 7 -l)(l + a)c-3-3( 7 -l)] , 



= 



HTs 



HT 



3a 



-a 



(7 - l)i 



it 

u v e 



e 



5 5- if 



(3.41) 



The eigenvalues that follow from Mnlcdh , using f|3.41[) and considering the high dissipative regime Q ^> 1, are still 
of the form obtained in the IS case, Eqs. (I3.33|) . (|3.34p and (I3.35|) . with coefficients as given by (|3.41[) . Obviously, one 
of the eigenvalue is still the same as the one obtained in the Eckart case, Eq. (|3.22l) . while the other two determine 
the stability condition for the NLCDH case, similar to Eq. (13.371) . 



1 



7 



-I - 



oo 1-39 



1 + 36/" 7- 1 1 + 36 2( 7 - 1) 1 + k 1 + 36 
For 7 = 4/3 the above equation gives 



1 + 36 



b > 



(3.42) 



1 



c + 4- 



3 00 1 — 36 



2 + 



b>0 



(3.43) 



1 + 36/ 1 + 36 2 1 + k 1 + 36 \ 1 + 36 

If we take the relaxation time as vanishing in Eq. (|3.43p . 6 = tH — > 0, we once again recover the result Eq. (|3.26[) 



IV. NUMERICAL RESULTS 



In this Section, we study numerically the system of equations for each of the three cases derived in the previous 
section. We will verify the corresponding stability conditions directly through the numerical time evolution of the 
corresponding dynamical systems. We will restrict to models of warm inflation that can be derived from quantum 
field theory. Moreover to develop the basic ideas, we will restrict to a dissipation coefficient of the form [§[ HH, HH 
T = CtfiT 3 /^ 2 , which is what is obtained from the two-stage mechanism in the low temperature regime. Although 
our analysis in this paper is restricted to this particular dissipative coefficient, it can easily be extended to other 
dissipative coefficients, with various other types studied in fi\, \l3HU§ . The bulk viscosity coefficient will have the form 
Cfc = CbT 3 , which is obtained for quantum fields in the high temperature regime [24| and also is the form generically 
considered in hydrodynamics.. This then corresponds to the case where c = 3 for the dissipative coefficient in Eq. 
(|3.3j) and I — 3 for the bulk viscosity, with and Cb (dimensionless) proportionality constants. Also for simplicity, 
we will analyze here the simplest case of a quadratic inflaton potential, 



with a constant relaxation time (g = g = g = g = 0in the constants defined in (|3.30|1 ). The extension to other types 
of potentials, such as a quartic potential or hybrid type potentials, does not offer additional difficulties and can be 
easily implemented. 
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In the example considered here, the stability conditions Eqs. 
NLCDH cases respectively, reduce to 



([3T26]) . (f3T38l) and (ET43)) . for the Eckart, IS and 



a 



Eckart 
stab 



3(l + cr)+4 + 3a-|- r ^- >0, 



2 1 



°stab 



r NLCDH 
°stab 



= 3 1 



K 



3 <7<7 



21 + K 

3 [ I 



+ 36 
2-39 
2 + 36 

17 



2 + 36 



86 



1 + 36 



2 

4- 



36 
3(T 



(1 + 



> , 

3 ad 1-36 
1 + 36 ~ 2 1 + k 1 + 36 



> 



(4.2) 

(4.3) 
(4.4) 



In all the numerical studies using the innaton potential Eq. (|4.1[) . we have kept fixed the values rti^ — v87rx 10 6 mp, 
the initial value for the dissipation factor Q = 100, the initial temperature T = 370 and 0(0) = 10.98 mp. The 
values for H(0) and </>(0) follow from the slow-roll conditions. These values correspond to a proportionality constant 
~ 1.61 x 10 8 for the dissipation term, which is a typical value found in the context of WI model building [40]. The 
value of the bulk viscosity coefficient C& is varied and also the value of the relaxation constant tH = 6, but observing 
that we are still in the region of validity of the thermal radiation bath at quasi-equilibrium, 6 < 1. 

By letting the system of equations evolve, we determine the critical value of Cb for which the stability conditions 
for each of the three theories studied here, Eqs. (|4.2I) . (I4.3[) and (|4.4I) . are violated. The corresponding results are 
given in Tab. |TJ 



e 


theory 


c b 




Eckart 


2232.94 


0.01 


IS 


2266.48 




NCLDH 


2300.04 




Eckart 


2232.94 


0.05 


IS 


2400.61 




NCLDH 


2568.42 




Eckart 


2232.94 


0.1 


IS 


2568.23 




NCLDH 


2903.79 



TABLE I: The critical values for the bulk viscosity constant Cb- The Eckart case is independent of ©, therefore the value for 
its critical Cb does not change. 



We note from the results of Tab. Q]that the values for the bulk viscosity constant Cb for which the stability conditions 
for IS and NCLDH cases are violated increases with respect to the Eckart case as 6 increases. In Tab. [IT] we give the 
corresponding differences in percentage. 



theory 


9 = 0.01 


e = 0.05 


e = o.io 


IS 


1.5% 


7.5% 


15.0% 


NLCDH 


3.0% 


15.0% 


30.0% 



TABLE H: The increase of the critical value of Cb for the causal theories with respect to the Eckart theory. 



Cnce we have the system evolving in time and also allowing the time dependence for the stability parameters 
^Jtat> art ' ^stab an d Citai? 131 ^ by starting from the initial conditions given above, we can explicitly check that the time 
where Eqs. (|4.2[) . (|4.3[) and (14.41) . are violated, is the time where both radiation energy density and the bulk pressure 
start to grow exponentially, as expected. In the Figs. [T] [2] and [3] we plot the stability parameters C , ^^f rt , Clf ah and 
C^atP DH along side those for the radiation energy density and the bulk pressure, as a function of time, for the case of 
6 = 0.01 and for the values of critical Cb shown in Tab. HI 

The results shown in Figs. [TJ [5] and [3] indicate that the time where the stability condition is violated corresponds 
to an inflection point in the radiation energy density and the bulk pressure evolutions. After that time, both the 
radiation energy density and the bulk pressure start to grow and soon after the dynamics become uncontrollable. The 
point where the stability conditions Eqs. (|4.2p . (|4.3I) and (|4.4I) . are violated, corresponds then to a turnover point in 
the evolution of the dynamical system of equations. 
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FIG. 1: The stability condition (a) and the results (b) for the radiation energy density, p r (black curves), and bulk pressure, 
II (green curves) , normalized by the total energy density. The dashed curves are for the Eckart case, the dotted curves are for 
IS and the dash-dotted curves are for NLCDH. In all cases Q = 0.01 and C b = 2232.94. 
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FIG. 2: The same as in Fig. [TJ but now for the critical value Cb = 2266.48. 



We can also notice from the results for the radiation and bulk pressure shown in Figs. [TJj, HJd and |3}3, that the 
causal theories always lead to a smaller radiation production when compared to the noncausal case. Among the two 
causal theories studied here, the NLCDH theory of Ref. [3l| gives a much smaller radiation production from bulk 
pressure effects than the IS theory. The differences between the causal theories of IS and NLCDH are also larger than 
the noncausal theory of Eckart as the relaxation time increases, which is clear from the results presented in Tabs. U 
and UU Even though the difference of the Eckart theory for the bulk pressure from that of IS is around the percent 
level for a relatively small relaxation time of tH = 0.01, it is twice that (in percentage) when the NLCDH theory 
is considered. This difference between the IS and NLCDH theories can easily be understood once we compare the 
Eqs. (|2.20[) and (|2.27l) and realizes that the relaxation time in the NLCDH theory appears with an additional factor 
two on the right-hand-side of the equation. Among the theories for the bulk pressure we have studied here, thus, the 
NLCDH case is the most robust in terms of stability. It allows for relatively larger bulk viscous pressures as compared 
to the Eckart and IS cases. 
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FIG. 3: The same as in Fig. [T] but now for the critical value Ct = 2300.04. 



V. BULK VISCOSITY COEFFICIENTS FROM QUANTUM FIELD THEORY 

The shear and bulk viscosities describe the properties of a system to return to equilibrium when displaced from it. 
As explained in [24j . at the level of particle physics processes, these viscosities are generally proportional to the mean 
free path, or equivalently time, of the relevant scattering process. The bulk viscosity is proportional to the mean free 
path for particle number changing processes in theories with breakdown of scale invariance. In contrast the shear 
viscosity is proportional to the two-body elastic scattering mean free path. Thus the bulk viscosity roughly has the 
form 



™4 

m t 



(5.1) 



where to is a characteristic measure of the violation from scale invariance in the theory and r is the mean free time 
between number changing inelastic scattering processes. 

Explicit expressions for the bulk viscosity have been calculated for a self-interacting \ y y A scalar field theoryjnodel 
in [24( for different temperature regimes. The obtained results for the bulk viscosity relevant for us here are [24| . 



5.5 x 10 
8.9 x l(T b 



4 "»Tj;i T) ln 2 [1.2465m y (T)/T] , m y < T < m y /X y 
A y V 3 ln 2 (0.064736A !; ), T^>m y /\ y , 



(5.2) 



where m y (T) is the scalar y field thermal mass, m 2 (T) = 



A y T 2 /24[l + 0(m y /T)], and 



m\ = 



ml{T) 



T 2 {dm 2 y {T)/dT 2 ) - to 2 - ^(A a )T 2 /48, where /3(X V ) = 3A 2 /(16tt 2 ) is the renormalization group /?-function. Note 
that even for a massless scalar field at tree-level, m y — and thus it is classically scale invariant, the scale invariance 
is broken quantum mechanically. In this case, the /3-function gives a measure of breaking of scale invariance. 

The characteristic relaxation time relevant for the bulk viscosity and that which also enters in the IS and NLCDH 
formulas for the bulk pressure, can be extracted from the result for the bulk viscosity given by Eq. (|5.2j) and the 
formal expression for it in the context of the Kubo formula for m y <C T and in the relaxation time approximation [4ll | , 



Cb - j! 



d 3 p t(uj p ] 
(2tt)3 W 2 



n{uj p )[l + n{uj p )} 



3 



2 2-2 



(5.3) 



where n(uj p ) = l/[exp(oj p /T) — 1] is the Bose-Einstein distribution, co p — yP 2 + m y {T) and v s is the speed of sound for 

the radiation bath of scalar y-particles. Using an on-shell approximation for the relaxation time, t(u p ) ~ r u> p /m y (T), 
where r = T(m y {T)) — constant, and the result for the speed of sound for a self-interacting scalar field in the high- 



temperature approximation m y <C T 24j, v 2 ~ 1/3 — 5to 2 /(127t 2 T 2 ), we obtain for the ratio C&/ T the result: 
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Cb ~ 4 1 f d 3 P 1 / 5p 2 1\ 

rhy T / 2T 

~ 187r 2 m w (T) n V m y( T ) 

Note that for the above results for the bulk viscosity to be applicable in warm inflation, it is required that the effective 
mass for the scalar y field be larger than the Hubble parameter, m y (T) 3> H. In this case curvature corrections to the 
quantum field expressions defining the bulk viscosity can be neglected and the Minkowski expression (15.21) applies. 
Also, as already pointed out in the previous sections, a quasi-equilibrium thermal radiation bath requires that the 
relevant relaxation time, set by r be also short compared to the Hubble time, t < 1/H. These conditions can be 
easily meet for warm inflation as verified in the next section. 




VI. MODEL BUILDING 



We have understood how to separate, in general, between the stable and unstable regimes. Moreover, we have 
studied the differences between the non-causal and causal descriptions of the bulk viscosity. In this Section, we are 
going to apply this knowledge to some generic supersymmetric (SUSY) models of warm inflation, namely, the chaotic, 
hybrid and hilltop (or new inflation) models. As an aside, recall that for single field, cold inflation models, cosmic 
microwave background radiation (CMB) data [42| is currently able to exclude many of the chaotic type of models, with 
potential V ~ <p n , such as with polynomial power n = 3 (at 95% CL) and higher powers (n > 3). Also typical hybrid 
inflation models are disfavored by the CMB data for having too large a spectral index n s > 1, while new inflation 
(hilltop) type cold inflation models seem to be more favored by the data [43} . However, these conclusions do not apply 
to warm inflation models [4(| , and in particular to large field inflation models with dissipation [26], [3t| , for which 
the predictions related to the primordial spectrum (spectral index, running of the spectral index, non-gaussianity 
parameter) arc different from those of standard cold inflation. In this paper we have only focused on the background 
evolution and so these observational details are not directly relevant. A detailed study of the primordial spectrum, 
that includes dissipation and viscous effects and the observational consequences, will be presented elsewhere. The 
goal of this Section is to study the parameter space of these basic inflation models in the presence of bulk viscosity. 

As discussed in previous sections, the bulk viscous pressure decreases the radiation pressure and so allows the source 
term creating radiation to be more effective. Thus assuming thermalization, it raises the temperature. Given that 
the dissipative coefficient depends on the temperature, the bulk viscosity enhances it; and therefore, the inflaton can 
slow- roll down its potential with lower values of the dissipative factor C^. Hence, we expect an enlargement of the 
parameter space in regions of low C^, where warm inflation is not allowed in the absence of bulk viscosity. In addition, 
in the regions of the parameter space where warm inflation is allowed without bulk viscosity, its effect is to produce 
more e-folds of inflation. 

One of the main issues that need to be faced in warm inflation are the quantum and thermal corrections to the 
inflaton potential, that tend to spoil the required flatness of the potential needed for inflation. We invoke SUSY to 
avoid the quantum corrections, which at the same time does not cancel the time non-local terms that give rise to the 
dissipative terms. The pattern of interactions leading to warm inflation is given by the superpotential: 



W — W(<&) + g$X 2 + hXY 2 , (6.1) 

where $, X and Y denote superfields. The regime most studied for warm inflation is where the X field, being coupled 
to the inflaton field, is massive during inflation, and therefore thermal corrections to the inflaton mass are Boltzmann 
suppressed as long as m x > T, where \ is the scalar component of the superfield X and m 2 x = 2g 2 <fi 2 . Even though 
the X field in (|6.ip is at low temperatures, we can still have the Y field at high temperatures, e.g., m y <C T, where 
y is the scalar component of Y and m y its mass. In this case, thermal dissipation of the inflaton field mediated by 
X decaying into Y is not Boltzmann suppressed [l]| and can lead to a warm inflation regime. We denote this regime 
as the low-T regime and it leads to a dissipative coefficient of the form [H, [l5| T = C^T 3 /(j) 2 . Other regimes 
and dissipative coefficients could also be studied, but we will develop the basic methodology in this paper with this 
example. 

Warm inflation imposes another condition on the temperature, T > H. This condition roughly separates warm 
inflation from cold inflation, since for T < H dissipation effects are negligible. We also must take into account the 
thermodynamical condition |n/p r | < 1, as the hydrodynamic descriptions of the bulk viscosity that we are using treat 
the viscous pressure as a perturbation to the equilibrium one. 
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Therefore, we define the available parameter space for warm inflation as the regions where the following conditions 
hold: 

1. eji = —H/H 2 < 1, which is the standard condition for the accelerated expansion, 

2. > p r , which prevents the radiation energy density to dominate, 

3. T/H > 1, which leads to non- negligible dissipation effects, 

4. 4>/T > 10, which is the low-T condition for g = 0(1), 

5. |II/p r | < 1, which is the condition for the hydrodynamic description to hold. 

These conditions need to hold for at least 40 e-folds to solve the flatness and horizon problems. In standard cold 
inflation, this value is taken to be 50-60 e-folds; nevertheless, the exact number depends on the scale of inflation and 
the initial post-inflation temperature. As warm inflation usually decreases both quantities, it is very reasonable to 
demand 40 e-folds. Condition 2 is controlled by the stability conditions found in the previous section. From system 
of equations (|3.12l) . we can relate the rest of the conditions with slow-roll parameters fl3.15p . Condition 1 during 
slow-roll simply is: 



e - = TTQ <1 - < 6 - 2 > 

The evolution of the ratio T/H, in the slow-roll regime, with respect to the number of e-folds is given by: 



dbx(T/H) 
H/pr is directly related to T/H, 



2(1 + a) /2 + 4Q \-Qrap 

-e - r\ + : - — V2e 



l + Q + 6Q(l + o-) V 1 



1 + 



(6.3) 



270C fc H 
7r 2 g» T 



(6.4) 



The evolution of <j>/T by: 



d\n(<j)/T) 
dN e 



-1 



[1 + Q + 6(5(1 + <r)] 



3 + 4q + (l + 2g)Q 
1 + Q 



2(l + q), + 3 + 2g+(5 + 4(T)Q ^V27 
1 + Q <j> 



(6.5) 



In addition to these equations, the slow-roll evolution of the field <f> is given by: 



dN e 1 + Q 

For completeness, we also show the evolution of Q: 



d<t>/m P _ V2i 



dQ 

dN P , 



Q 



l + Q + QQ{l + a) 



10 1 



6(l + (r)ij + 8 1 + 



3 \ m P 
-a ) —V2e 



(6.7) 



These results generalize the ones obtained in [40j for the case with no bulk viscosity. Note the difference in the 
notation between our a and the defined there, which have replaced here by V2e(cf>/mp). Nevertheless, as we have 
shown in (|3 . 24[) . |<r| < 1 and then the results in concerning whether the conditions increase or decrease during 
the evolution are still valid. 

The last question before entering in the particular details of each model is how we are treating the bulk viscosity. 
In this Section we will use the non-causal description of the bulk viscosity, i.e. the Eckart theory, and place limits on 
the validity of this approximation. As discussed previously, the Eckart description is a good approximation for low 
values of O = tH. The Hubble parameter is given by (|2.7I) and the relaxation time r is obtained from Eq. (15.41) . 
Using for example the first expression for the bulk viscosity in Eq. (|5.2p (the second expression for the bulk viscosity 
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in Eq. (|5.2[) can be easily seen as obtained from the first, when neglecting the temperature independent terms in 
m y (T) and in rh y ), we obtain that 



T « 9 .77 x 10 6 ^/, (6.8) 

where for the superpotential (|6.ip . we have that Aj, = 6h 2 . Based on the previous discussions, we consider the 
Eckart approximation to be valid when < 10~ 2 , which translates into (using that in the high-temperature limit 
m y (T) w hT/2) 



h > 10 - . (6.9) 

One of the conditions for warm inflation is that T/H > 1, therefore, we can easily arrange the condition (I6.9P to be 
satisfied when deep in the warm inflationary regime, particularly in the strong dissipative regime, which can also allow 
for perturbative values for the coupling h. The effect of including a causal description of the viscosity is to produce a 
lower value of the bulk viscous pressure than in the non-causal case with the same coefficient Cb- The consequence is 
a shift around the 5% level in the entire parameter spaces shown in the next subsections to higher values of CV This 
is explicitly verified below for the specific inflaton potential models we have studied. 



A. Chaotic model 



We first consider a chaotic inflation potential of the form: 



V(<f>) = (6.10) 

where we have used A = 10~ 14 . However this parameter is only relevant for the amplitude of the power spectrum, 
which we are not interested in here. Using again that T = C^T 3 /^ 2 , for the potential (|6.10p the slow-roll parameters, 
Eq. p,15p . are given by 



"= 12 (T) 2 ' «-§«■ "-§«■ <fm) 

Therefore, the value of the field decreases during inflation, meanwhile, the dissipative ratio Q and T/H both increase. 
Hence, once the condition on T > H is fulfilled initially, it is always satisfied. The ratio (j>/T decreases, but we have 
checked that it always remains above 10 as long as the other conditions are fulfilled. The parameter en increases 
during inflation and, as a consequence, warm inflation ends when en = 1- Finally, the condition > p r is controlled 
by the stability condition (|3.26l) . It is only necessary to check that the stability condition is positive at the beginning 
of inflation, as it does not change sign during the evolution. This last statement is true for all the models studied. 

The available parameter space is shown in Fig. For completeness we have included the parameter space excluded 
for different values of the hydrodynamic condition |n/p r |, namely 0.1, 0.5 and 1. We observe that the enlargement of 
the parameter space in regions of low C$ is not very efficient. In particular, the minimum value of C^ is reduced from 
2.1 x 10 6 up to 1.6 x 10 6 . These values have been confirmed by using the NLCDH description of the bulk viscosity, 
Eq. (|2.27p . with a constant r fixed by imposing initially the values = 0.01, 0.9. The NLCDH description reduces 
the bulk viscous pressure associated to a Cb value when grows, which means that higher values of Cb are allowed 
before the condition |n/p r | is violated. However, at the same time, for the same value of C^, higher values of Cb are 
required to avoid the T/H < 1 exclusion region, therefore the effects cancel each other and the minimum value of 
is independent of the initial value 0. 

We have also found that the condition |II/p r | is the most restrictive one in almost the full parameter space and 
that the instability regime studied in the previous sections is far beyond the limit imposed by this condition. In the 
analysis we have fixed the initial values such that the exclusion regions are the least stringent, that is, the upper 
region is as high as possible and the bottom one, as low as possible. We have fixed initial conditions in this way for 
the three models studied. 
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In addition, we separate with black lines the regions where the dissipative ratio at horizon crossing Q* is always 
greater than one from the regions where it is always less than one. In the region between them, Q* can be either 
greater or less than one, depending on the initial value of 4>. In the regions that were not allowed in the absence of 
bulk viscosity, its main role is to produce enough e-folds of inflation. In regions allowed with no bulk viscosity, the 
total number of e-folds is increased. We quantify this effect in terms of the percentage difference AiV e , defined as 



^ybulk ^yno-bulk 

AiV e = 6 Ar . e .. x 100, (6.12) 

ATno-bulk v ' 

where ulk is the maximum number of e-folds obtained with bulk viscosity for a certain C$ and N^°~ bulk is the 
equivalent without bulk viscosity. The results are shown in Fig. [5] 




FIG. 4: Parameter space for the chaotic model. The green regions are excluded because of the violation of the condition written 
in the plot. The lines separate regions where Q* < 1 and Q* > 1 respectively. In the region between them, we can have both 
Q, < 1 and Q, > 1 for different values of 0(0) . 

The bulk viscosity can enhance the number of e-folds through two mechanisms. The first one is, for a given initial 
value of the field, to reduce the redshift of the radiation energy density by decreasing its total pressure. This effect 
produces an increase in the value of Q, which goes as f?J^ and is related to the number of e-folds via 



w .= r' 3g2( ' +(a #, (6.i3) 

where 0(0) is the initial value of the field and (j) e nd, the value of the field at the end of inflation. Hence, the increase in 
Q leads to a larger number of e-folds. However, we have checked that this mechanism is subdominant in the quartic 
potential, providing an efficiency up to 2%. 

The second mechanism allows to increase the initial value of the field. From Eq. (|6 . 1 3[) it can be seen that this 
produces more e-folds by increasing the integration interval. In the absence of bulk viscosity, there is an upper limit 
on the value of the field because of the condition T/H > 1. As the bulk viscosity increases the value of T, it pushes 
upwards this limit and, hence, it is possible to use larger values of the field. However, there is a bound to this 
effect, provided by the condition |n/p r | < 1, which translates into a new upper limit to </>(0). The decrease in this 
effect shown in Fig. [5] is because the limits on the initial value of the field grow at different rates with C^, with the 
former growing faster, which causes the difference on the limits to decrease with and this, in turn, reduces the 
difference between the viscous and non- viscous cases. From Eq. (|6.4I) we can obtain the value of Cb that maximizes 
this mechanism. Using the values |n/p r | = 1, T/H = 1 and <?* = 225.78 as an example, we find that Cb = 8.25. This 
argument is model independent, so that we find the same value of Cb in the three models studied and independently 
of the value of C$. 



B. Hybrid models 



We consider now small field models of inflation with an inflationary potential given by: 
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FIG. 5: Maximum enhancement in the number of e-folds for different values of Cf, for the model with quartic potential. The 
maximum value is obtained for Ct=8.25. 



V = V 



1+1 

2 \ rap 



(6.14) 



where Vq is the scale and 5 a model parameter. Here we have used Vq = 10~ s m P and study the dynamics for two 
representative values for the parameter 8. The slow-roll parameters, in the case of the inflaton potential given by 
Eq. (|6.14l) . are now given by 



2 \mp 



= -2n. 



(6.15) 



During the evolution the value of <j> decreases, while the value of Q increases. The evolution of the value of T/H 
depends on the value of Q: for Q < 1, T/H increases and for Q > 1, it decreases; cjy/T always decreases, and e# is 
not relevant in this model because it is suppressed by a factor (0/mp) 2 , which is usually very small due to the fact 
that <j) -C m p . Hence, e# is always below one. As a result, inflation ends because the conditions imposed either on 
T/H or (j>/T are violated, or because the field reaches its critical value. 

In Fig. El the parameter space of the hybrid models with S = 0.1 (left panel) and 8 = 10 (right panel) are plotted. 
In the left panel, as in the previous model, the black lines separate regions with different value of the dissipative ratio 
at horizon crossing. However, in the right panel, the dissipative ratio is always above one. For a higher value of 8 we 
can maintain tj/(1 + Q) below one only for > 1. In the left panel the bottom excluded region is forbidden for the 
same reason as in the chaotic model. Nevertheless, in the right panel we find that the excluded region is forbidden 
by the <j>/T > 10 condition. This is caused again by the higher value of rj. The parameter 8 in Eq. (|6. 14[) measures 
the curvature of the potential. Thus, for higher values of 8 the field evolves faster. As a result, the condition on <j)/T 
is reached first than in the case for smaller values for the parameter 8. 

The minimum value of in this case of inflation with the hybrid type of potential, Eq. (I6.14[) . is reduced from 
3.5 x 10 4 up to 2.6 x 10 4 , for 8 — 0.1, while for 5 = 10 it is reduced from 5 x 10 4 up to 4.1 x 10 4 . Making use of 
the NLCDH description of the bulk with a constant r fixed by imposing the initial values 6 = 0.001, 0.9, we found 
that the minimum value is independent of the initial choice of G. The effect on the number of e-folds is shown in 
Fig. [7] As in the previous case, the constant field mechanism is subdominant, with an efficiency of around a 3% for 
the 8 = 10 case and a 6% efficiency for the 8 — 0.1 case. 

C. Hilltop models 

We now consider new inflation hilltop type of models, which are characterized by an inflaton potential given by 



V = Vn 



2 \mp 



(6.16) 
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FIG. 6: Parameter space for the hybrid models with 5 — 0.1, 10 The green regions are excluded because of the violation of the 
condition written in the plot. In the left plot the lines separate regions where Q* < 1 and Q* > 1 respectively. In the region 
between them, we can have both Q* < 1 and Q, > 1 for different values of 4>{0). In the right plot, Q, is always above 1. 
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FIG. 7: Maximum enhancement in the number of e- folds for different values of Cf, in the case of the inflaton potential given 
by Eq. (|6.14|) . The maximum value is obtained for Cf, = 8.25. 



where the dots account for higher-order terms and Vq = 10~ s rrip. This is a potential similar to the two previous ones, 
but with a negative squared inflaton mass. The slow-roll parameters are still given by those in Eq. (|6.15l) . but with 
the change 6 — > —5. In these models <fi, <f>/T, T/H and ch increase during the evolution, while Q decreases. Inflation 
ends when the field reaches a large enough value, so that higher-order terms in the potential start contributing and 
ch becomes greater than one. 

The parameter space for the hilltop model Eq. (|6 . is shown in Fig. [SJ for the cases of 6 = 0.1 (left panel) and 
5 = 1 (right panel). Once again, in the left panel the black curves separate regions with Q* greater or less than one 
at horizon crossing, and in the right panel, the dissipative ratio at horizon crossing is always greater than one due to 
the large value of r\. 

Now, the minimum value of is reduced from 3.3 x 10 up to 2.5 x 10 4 for 8 = 0.1 and from 5.9 x 10 up to 
5.1 x 10 4 for 6 = 1. We have checked these values with the NLCDH description of the bulk with a constant r fixed 
by imposing the initial values 9 = 0.001, 0.9. The effect on the number of e-folds is shown in Fig. [9] As in the other- 
cases studied, the constant field mechanism is subdominant with efficiencies around a 5% and a 3% for the <5 = 1 and 
the 5 = 0.1 cases respectively. The initial value mechanism works reversely compared with the previous potentials. 
Here it allows to use lower initial values of the field, however, as the value of the field grows in this case, this reduction 
implies an increase of the integration interval in Eq. (|6.13[) . In addition, note that in this potential, the lower the 
value of the field, the larger the value of H and thus, the lower is the ratio T/H . This argument also applies to the 
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|II/p r | condition, therefore, there are lower limits to </>(0) rather than upper ones. 




FIG. 8: Parameter space for the hilltop models with S = 0.1, 1. The green regions are excluded because of the violation of the 
condition written in the plot. In the left plot the lines separate regions where Q* < 1 and Q* > 1 respectively. In the region 
between them, we can have both Q, < 1 and Q» > 1 for different values of <j)(0). In the right plot, Q* is always above 1. 
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FIG. 9: Maximum enhancement in the number of e-folds for different values of C^, for the case of the inflaton potential 
Eq. (|6.16[) . The maximum value is obtained for Cb = 8.25. 



D. The bulk viscosity coefficient 



Let us now briefly discuss the expected values for the bulk viscosity coefficient for a warm inflation model described 
by Eq. (I6.1|) . In the context of the model described by Eq. (|6.1|) . for the relevant temperature regime for warm 
inflation, characterized by m x > T and m y <C T, the dominant contribution to the bulk viscosity comes from the 
radiation thermal bath composed of the high-temperature particles, i.e. the y particles and given by Eq. (|5.2|) . We 
can see from the expression for the bulk viscosity Eq. (|5.2j) that the larger values for the bulk coefficient Cb = (b/T 3 
appears in the intermediate temperature regime, m y <T< m y /X y , where the y particles are already in the high- 
temperature regime, but the temperature is still not too high, such that in the thermal mass m y (T) the temperature 
corrections are subdominant. In this case, neglecting the thermal corrections to the mass, we get the estimate for Cb, 



1 ml 



C b w 5.5 x lO^^ln 2 (1.2465*4/0 , (6.17) 
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recalling that for the model (|6.ip . X y — 6ft. 2 . Taking m y /T ~ 0.1, we obtain 70 < C*6 < 1.8 x 10 4 , for values of h 
between 0.1 and 0.2. This is in the absence of further decay modes, which would increase even more the estimates 
for Cb (the bulk viscosity coefficient is proportional to the radiation bath field degeneracy) . These values are already 
within the window of values of the viscosity coefficient observed by the results in Figs. HI |5]and[Sl which allows warm 
inflation with smaller dissipation as a consequence of including a bulk viscous pressure. 

VII. CONCLUSIONS 

In this work we have examined the stability of the dynamics of warm inflation when bulk viscous effects are included. 
The noncausal theory for the bulk viscous pressure, given by the Eckart hydrodynamics theory, was studied along 
with two other causal theories, the IS linear theory and the NLCDH theory, which is a simpler causal theory for the 
bulk viscosity, yet nonlinear. The stability of the resulting dynamical system for warm inflation has been studied 
before by the authors of Refs. [2(|, in the absence of bulk viscous pressure and [26|, by including the effect of a viscous 
radiation bath for only the Eckart case. In this work we have extended these previous works by also carrying out the 
analysis of stability for the causal theories. We have also explored model building in the context of warm inflation 
with bulk viscous pressure. 

By including the analysis of stability for the causal theories of IS and NLCDH, it becomes possible to clarify on 
the differences between these two theories and with that of the noncausal Eckart theory. We have seen significant 
differences for the radiation production in each of these different theories as the relaxation time of the radiation 
fluid increases. Among the three theories for the bulk pressure we have studied, the NLCDH case proved to be the 
most robust of them as far as stability is concerned. We have shown that it allows for relatively larger bulk viscous 
pressures as compared to the Eckart and IS cases. The noncausal theory is the one that most quickly can develop 
instabilities. Our results also suggest the significance and differences there might be in other nonlinear generalizations 
of theories for the viscous hydrodynamics of radiation fluids in cosmology, when relatively larger relaxation times are 
involved, but still within the range of validity of thermal equilibrium for the radiation bath, < 1. Recalling that 
the Eckart and the IS theories have a long history of use in cosmology in general [22], [23| , other nonlinear realizations 
of viscous hydrodynamics, in particular for the bulk pressure, may lead to significant differences in results, similar to 
the comparisons we found here for our problem. 

In regards model building, we have shown that accounting for bulk viscous pressure effects in the radiation fluid can 
relax the requirements on the magnitude of the dissipation coefficient for the inflaton field, especially for a large bulk 
viscosity coefficient. This range of bulk viscosity coefficients can be realized within the regime of stability requirements 
in warm inflation, and this range is within reach of realistic model parameters. The upshot is, accounting for bulk 
viscosity effects makes warm inflation more robust. In particular warm inflation is more readily realized and with less 
constraints when bulk viscosity effects are accounted for versus when they are ignored. Although our examination of 
warm inflation was only for one particular form of the dissipative coefficient, we believe this conclusion will hold in 
general. Our results point to the importance of considering viscous effects in future studies of warm inflation. We 
also believe our results can have more general implications for systems involving viscous radiation fluids, in particular 
for cosmology. Very little has been discussed in the literature on the significance to cosmology of actual expressions 
of bulk viscosity computed from underlying particle physics models. This is one of the few papers that has examined 
this question for the specific case of warm inflation dynamics. However the methodology developed here can easily be 
extended to examine the same sorts of bulk viscosity effects in other relevant cosmological epochs, such as during the 
reheating phase after inflation in cold inflation dynamics and during the radiation and matter dominated expansion 
regimes. 

Acknowledgments 

M.B-G. and R.C. are partially supported by MICINN (FIS2010-17395) and "Junta de Andalucia" (FQM101) A.B. 
acknowledges support from the STFC. R.O.R is partially supported by research grants from Conselho National de 
Desenvolvimento Cienti'fico e Tecnologico (CNPq) and Fundacao Carlos Chagas Filho de Amparo a Pesquisa do Estado 
do Rio de Janeiro (FAPERJ). CSV. is supported by Coordenagao de Aperfeigoamento de Pessoal de Ni'vel Superior 
(CAPES). R.O.R. and CSV. would also like to thank T. Koide for discussions concerning the hydrodynamical 



27 



equations with bulk viscosity. 



[1] A. H. Guth, Phys. Rev. D 23 (1981) 347; K. Sato, Phys. Lett. B 99 (1981) 66; A. Albrecht and P. J. Steinhardt, Phys. 

Rev. Lett. 48 (1982) 1220; A. D. Linde, Phys. Lett. B 108 (1982) 389. 
[2] A. D. Linde, Phys. Lett. B 129 (1983) 177. 

[3] A. Berera and L. Z. Fang, Phys. Rev. Lett. 74 (1995) 1912 [arXiv:astro-ph/9501024]. 

[4] A. Berera, Phys. Rev. Lett. 75 (1995) 3218; [arXiv:astro-ph/9509049]. Phys. Rev. D 54 (1996) 2519; [arXiv:hep- 

th/9601134]. Phys. Rev. D 55 (1997) 3346. [arXiv:hep-ph/9612239]. 
[5] L. Z. Fang, Phys. Lett. B 95 (1980) 154 ; I. G. Moss, Phys. Lett. B 154 (1985) 12; J. Yokoyama and K. I. Maeda, Phys. 

Lett. B 207 (1988) 31. 

[6] E. A. Calzetta and B.-L. B. Hu, Nonequilibnum Quantum Field Theory (Cambridge University Press, Cambridge, 2008). 
[7] A. Berera, M. Gleiser, R. O. Ramos, Phys. Rev. D 58 (1998) 123508. [arXiv:hep-ph/9803394 [hep-ph]]. 
[8] A. Berera, I. G. Moss and R. O. Ramos, Rept. Prog. Phys. 72 (2009) 026901. [arXiv:0808.1855 [hep-ph]]. 
[9] D. Boyanovsky, D. Cormier, H. J. de Vega and R. Holman, Phys. Rev. D 55 (1997) 3373 [hep-ph/9610396]. D. Boyanovsky, 
D. Cormier, H. J. de Vega, R. Holman and S. P. Kumar, Phys. Rev. D 57 (1998) 2166 [hep-ph/9709232]. 
[10] M. Gleiser and R. O. Ramos, Phys. Rev. D 50 (1994) 2441 [hep-ph/9311278]. 

[11] N. Barnaby, E. Pajer and M. Peloso, Phys. Rev. D 85 (2012) 023525, [arXiv:1110.3327 [astro-ph.CO]]; D. Lopez Nacir, 

R. A. Porto, L. Senatore and M. Zaldarriaga, JHEP 1201 (2012) 075, [arXiv:1109.4192 [hep-th]]. 
[12] A. Berera, M. Gleiser, R. O. Ramos, Phys. Rev. Lett. 83 (1999) 264. [arXiv:hep-ph/9809583] . 
[13] M. Bastero-Gil, A. Berera and R. O. Ramos, JCAP 1109 (2011) 033. [arXiv:1008.1929 [hep-ph]]. 
[14] M. Bastero-Gil, A. Berera, R. O. Ramos and J. G. Rosa, [arXiv: 1207. 0445 [hep-ph]]. 
[15] I. G. Moss and C. Xiong, [arXiv:hep-ph/0603266]. 

[16] A. Berera, R. O. Ramos, Phys. Rev. D 63 (2001) 103509. [arXiv:hep-ph/0101049]. 
[17] A. Berera and R. O. Ramos, Phys. Lett. B 567 (2003) 294 [arXiv:hep-ph/0210301]. 

[18] A. Berera and R. O. Ramos, Phys. Lett. B 607 (2005) 1 [arXiv:hep-ph/0308211]; Phys. Rev. D 71 (2005) 023513 [arXiv:hep- 
ph/0406339]. 

[19] L. M. H. Hall and I. G. Moss, Phys. Rev. D 71 (2005) 023514 [arXiv:hep-ph/0408323]. 
[20] I. G. Moss and C. Xiong, JCAP 0811 (2008) 023 [arXiv:0808.0261 [astro-ph]]. 
[21] S. Weinberg, Gravitation and Cosmology, (New York, NY, Wiley, 1972). 

[22] G. L. Murphy, Phys. Rev. D 8 (1973) 4231; Y. B. Zeldovich, Pisma Zh. Eksp. Teor. Fiz. 12 (1970) 443; B. L. Hu, Phys. 

Lett. A 90 (1982) 375; I. Waga, R. C. Falcao and R. Chanda, Phys. Rev. D 33 (1986) 1839; R. Maartens, Class. Quant. 

Grav. 12 (1995) 1455; W. Zimdahl, D. Pavon and J. Triginer, Helv. Phys. Acta 69 (1996) 225 [arXiv:astro-ph/9608060]. 

W. Zimdahl and D. Pavon, Phys. Rev. D 61 (2000) 108301; M. Giovannini, Class. Quant. Grav. 22 (2005) 5243 [arXiv:astro- 

ph/0504655]; S. del Campo, R. Herrera, D. Pavon, Phys. Rev. D 75 (2007) 083518 [arXiv:astro-ph/0703604]. 
[23] J. -S. Gagnon and J. Lesgourgues, JCAP 1109 (2011) 026 [arXiv:1107.1503 [astro-ph.CO]]; I. Brevik, E. Elizalde, S. Nojiri 

and S. D. Odintsov, Phys. Rev. D 84 (2011) 103508 [arXiv: 1107.4642 [hep-th]]; O. F. Piattella, J. C. Fabris and W. Zimdahl, 

JCAP 1105 (2011) 029 [arXiv:1103.1328 [astro-ph.CO]]; C. Bogdanos, A. Dimitriadis and K. Tamvakis, Phys. Rev. D 75 

(2007) 087303 [arXiv:hep-th/0611094]. 
[24] S. Jeon, Phys. Rev. D 52 (1995) 3591 [arXiv:hep-ph/9409250] ; S. Jeon and L. G. Yaffe, Phys. Rev. D 53 (1996) 5799 

[arXiv:hep-ph/9512263]. 
[25] H. P. de Oliveira and R. O. Ramos, Phys. Rev. D 57 (1998) 741 [arXiv:gr-qc/9710093]. 

[26] S. del Campo, R. Herrera, D. Pavon and J. R. Villanueva, JCAP 1008 (2010) 002 [arXiv:1007.0103 [astro-ph.CO]]. 

[27] J.P. Mimoso, A. Nunes and D. Pavon Phys. Rev. D 73 (2006) 023502 [arXiv:gr-qc/0512057] 

[28] C. Eckart, Phys. Rev. 58 (1940) 919. 

[29] W. Israel, Ann. Phys. (NY.) 100 (1976) 310. 

[30] W. Israel, J. Stewart, Ann. Phys. (NY.) 118 (1979) 341. 

[31] G. S. Denicol, T. Kodama, T. Koide and P. .Mota, J. Phys. G 36 (2009) 035103. [arXiv:0808.3170 [hep-ph]] 
[32] R. Maartens and V. Mendez, Phys. Rev. D 55 (1997) 1937 [arXiv:astro-ph/9611205]. 
[33] P. Ilg and H. C. Ottinger, Phys. Rev. D 61 (2000) 023510. 

[34] M. Bastero-Gil, A. Berera and R. O. Ramos, JCAP 1107 (2011) 030 [arXiv: 1106.0701 [astro-ph.CO]]. 

[35] R. Maartens, [arXiv:astro-ph/9609119]; R. Maartens, J. Triginer, Phys. Rev. D56 (1997) 4640 [arXiv:gr-qc/9707018]. 

[36] W. A. Hiscock and L. Lindblom, Annals Phys. 151 (1983) 466. 

[37] L. Landau, E.M. Lifshitz, Fluid mechanics (Addison Wesley, Reading, MA, 1958). 

[38] A. R. Liddle, P. Parsons and J. D. Barrow, Phys. Rev. D 50 (1994) 7222 [arXiv:astro-ph/9408015]. 

[39] R. O. Ramos and L. A. Silva, work in progress. 

[40] M. Bastero-Gil, A. Berera, Int. J. Mod. Phys. A 24 (2009) 2207 [arXiv:0902.0521 [hep-ph]]. 
[41] P. Chakraborty and J. I. Kapusta, Phys. Rev. C 83 (2011) 014906 [arXiv:1006.0257 [nucl-th]]. 
[42] E. Komatsu et al, Astrophys. J. Suppl. 192 (2011) 18 [arXiv:1001.4538 [astro-ph.CO]]. 

[43] D. Boyanovsky, C. Destri, H. J. De Vega and N. G. Sanchez, Int. J. Mod. Phys. A 24 (2009) 3669 [arXiv:0901.0549 
[astro-ph.CO]]. 



